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Abstract. This paper studies a generalization of the Curie- Weiss model (the Ising model 
■^^i on a complete graph) to quantum mechanics. Using a natural probabilistic representation 

of this model, we give a complete picture of the phase diagram of the model in the pa- 
rameters of inverse temperature and transverse field strength. Further analysis computes 
the critical exponent for the decay of the order parameter in the approach to the critical 
curve and gives useful stability properties of a variational problem associated with the 

^ 1. Introduction 

^ S ^ One of the simplest classical systems exhibiting phase transition is the Curie- Weiss model. 

In this model, N Ising spins (i = {cri = ±1; / = !,..., N}, interact via the Hamiltonian 

1 ^ 

^ HN{^) = -j^L'^''^j' (1-1) 

O 

^ where the normalization by ^/n makes Hj^ a quantity of order N. (Generalizations to 

multiple-spin interactions may also be considered; cf Sect. |3] However, for the time 



^ being, will suffice.) 

00 As is well known for the spins distributed according to the measure cx 

O Q-^HN{a)^ j jj-^ ^Yxe canonical ensemble, as N — > oo the law of the empirical mean 

_ ^ m^[a) = N^^ J^i ai converges to a mixture of point masses at rk.m^ where m^ — 

^ is the so called spontaneous magnetization. The phase transition in this model is mani- 

^ festly seen from the observation that m*(j6) = for j6 < /5c while m*(j6) > for /5 > fic- 
The function /3 i— > is in fact the maximal non-negative solution of the equation 

m = tanh(/3nj), (1.2) 

while /3c = 1 is the so called critical temperature. 
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The main reason why the Curie-Weiss model is so approachable is the fact that the 
Hamiltonian is, to within an additive constant, equal to —2Nmf^{a)^. This permits 
a very explicit expression for the law of whose concentration properties are then 
readily controlled by straightforward large-deviation arguments. However, this sim- 
ple strategy breaks down once inhomogeneous terms (i.e., those not invariant under 
exchanges of the spins) are added to the Hamiltonian. One example where this hap- 
pens is the Curie-Weiss system in random external field where the term hjCTi, with hi 
sampled from an i.i.d. law with zero mean, is added to Hjv. While rigorous analysis is 
still possible in this case, the technical difficulties involved are more substantial. Signifi- 
cantly more complex is the Sherrington-Kirkpatrick version of where the term (7,c7"y 
is weighed by a (fixed) random number }x,ij that has been sampled from a symmetric 
distribution on R. This model possesses a beautiful underlying structure 1201 [23 1 which 
has been harnessed mathematically only very recently [17, 25. 26]. 

The goal of this paper is to study another natural generalization of the Curie-Weiss 
model, namely to the realm of quantum mechanics. Here the classical spin variables CTj are 
understood as eigenvalues of z-component of the triplet of Pauli matrices {d'^-^\ a^y\ a^^"*)- 
the generators of su(2) — acting on the one-particle Hilbert space J^i = SPAN{ | +) / 1 — ) }• 
The configuration space is replaced by the product space = <^fLi '^i- Classical Ising 
cofiguartions a G {±1}^ generate an orhtonormal basis \gi) = ^ \cri) of J^n- If ^ is an 
operator on (e.g. A = a^^^ or A = a^^^) then its copy A/ acts on i-th spin component 

of a product vector cp) = '^^^i \(pi) G J^n via 



Ai = \(pi)®---®A \(pi) (g) ■ ■ ■ (g) |(Pn) . (1.3) 

There are at least two natural ways to introduce quantum effects into the Curie-Weiss 
model. Either one may make the interaction term isotropic — this corresponds to the 
quantum Heisenberg model — or one may consider an external transverse field. Here 
we focus on the latter situation: The Hamiltonian is now an operator on .^jq defined by 



-1 N N 

= E - Zihaf + A^W). (1.4) 

i,j=l i=l 



Instead of Gibbs probability measures one studies a KMS state { — )^f which is a pos- 
itive linear functional on the C*-algebra of operators on J^j-j defined by 

M./3_ Tr(Ae-/'^^) 

As before, the parameter j5 plays the role of inverse temperature while A, which corre- 
sponds to the strength of an external field, determines the overall strength of the quan- 
tum perturbation. 
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In particular, ( |1.5| l gives rise to a probability distribution ji^f on classical spin con- 
figuartions (7 G {±1}^, 

?^a; [^] = Tr{e-f^nly (1-6) 

and one can try to read a signature of phase transition in terms of concentration proper- 
ties of the latter as N tends to infinity. 

Despite this relatively clean formulation, the quantum nature makes this model very 
different from the classical one. Indeed, the Hamiltonian is not simultaneously diago- 

(z) 

nalizable with the &■ 's for A 7^ 0, so the very notion of "value of spin at /" is apparently 
lost. In fact, as we explain below, a path integral approach or, equivalently, a stochas- 
tic geometric representation of ji^f reveals that relevant large deviation and measure 
concentration analysis in the genuine quantum case A > should be lifted to infinite 
dimensions. 

We refer to (Si |6l HI HU |22| where the ideas of stochastic geometric representation in 
question were originally developed, as well as to a recent review fT8\. The bottom line 
is the representation ji^f- Let be the circle of circumference fi. There exists a shift 

invariant probability measure ji^ on piece-wise constant trajectories a : Sp {il}/ 
such that 



N 

! = 1 



N Jg^llmf,{tf+hmN{t)]dt J 



(1.7) 



wherea{t) = {cri{t), ... ,o-fq{t)) and m^{t) = N"^ cr,(i). 

With such a representation in place, we then apply methods of large-deviation the- 
ory to derive the leading-order N — > 00 asymptotic of these expectations. As can be 
expected, many physical properties of the quantum system may be gleaned from the 
properties of the minimizer of the corresponding variational problem. 

The variational problem (or rather dual thereof) turns out to have intrinsic features 
which allow us bring methods of FX percolation to bear. Many quantitative charac- 
teristics of the system will be determined explicitly (at least in the limit N — > 00). In 
particular, we obtain full control of the phase diagram and the stability near optimiz- 
ers for variational problem. Moreover we give an explicit characterization of the critical 
exponent for the decay of the optimizer as one approaches the critical curve. 

Let us mention existing work on issues encompassing aspects of the present paper. 
The paper |10 | addresses general mean field quantum spin systems in their C* algebraic 
representation. Relying on an operator theoretic version of deFinetti's Theorem (i.e. 
Stormer's Theorem), they derive a 'mean field equation' for the extremal states of the 
system and formulate variational problem which these states must solve. More recently, 
the preprint [9J uses a slightly different path integral representation from ours to derive 
a qualitative probabilistic description of the solutions to these 'mean field equations'. 
However, neither work attempts a detailed analysis of the Quantum Curie-Weiss Model. 
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The remainder of this paper is organized as follows. In Section|2]we review the prob- 
abilistic formulation of our problem and use this opportunity to set notation. In Section 
|3] we formulate the main results of this paper. Subsequent sections are devoted to the 
proof of these results. 



2. Stochastic Geometry of the Mean Field Transverse Ising Model 



Let us explicitely describe the one-circle measure f/^ in ( |1.7| |. It should be noted that a 
natural approach, exploited in the papers Aizenman-Klein-Newman [5J and Campanino- 
Klein-Perez ||8l among others is to realize the Gibbs state f/^ = (■)^ which we construct 
directly here, as the strong coupling limit of a sequence discrete Ising systems on Z/ nZ. 
As a result of this approximation we get a weak form of ferromagnetic correlation in- 
equalities for (■)^ ^ for free. We well rely on this below without further mention. 

We begin by introducing a convenient notation for expectation values. Given a prob- 
ability measure P on a sample space Q and an integrable function / : Q — > R, let us 
denote the expectation value of / by P (/) . 

Let P^ be the distribution of the Poisson point process ^ C S^ of marks(which we 
visualize as puntucres) on the circle S^ with arrival intensity A. Given a realization of ^ 
let us say that a classical piece-wise constant trajectory c7 : S^ i— > {±1} is compatible with 
^; cr ~ ^, if jumps of cr{-) occur only at arrival times of Note that we do not require 
that cr changes sign at each arrival of ^. Consider now a joint probability distribution. 



0f jd?,d(7] cxP$[d^]e''/'-(Odf^^^^^^. 



(2.1) 



Here <^{£r~^} gives mass one to a if it is compatible with ^ and gives zero mass otherwise. 
In the sequel we shall supress the sub-index h whenever h = 0. Our one-circle measure 
fi^ is just the t7"-marginal of O^. Clearly, for h ^ 0, cr-marginals ji^j^ of can be 

recovered from /^^ as exponential tilts by a constant magnetic field h, 



FA,h l^^] V? 



Summing with respect to compatible c-trajectories in (2.1 1 we recover the ^-marginal 
Q^^ of O^: For each realization of ^ the punctured circle S^ \ ^ is a disjoint union of 
finite number #(^) of connected components. 



7=1 



Then, 



Ql. m - n 



7=1 
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Note that in the case h = 0, 

In general, ;,| is still a family of Fortuin-Kasteleyn type random cluster measures. 
In particular, they possess two specific properties which we would like to stress: 

1) Stochastic domination. For each h there exists an intensity rj = r]{h,\) of the Poisson 
process of arrivals ^ on S^, such that 

^ P? (2.3) 

We refer to ||3l|5l for a startegy of proving statements like (23 1, which is based on a 
comparison with independent Bernoulli percolation via the strong FKG inequality. 

2) Edwards-Sokal representation. The spin marginal /^^^ could be recovered from Q^^ 
by the following two step procedure: First sample f from ^, and then independently 
paint each cormected component 7 of \ ^ into ±1 with probabilities e''l^l / (e'^'^l + e^'^'^l ) 
and e-''l^l/(e''l^l + e^'^l^l) respectively 

In the sequel we shall repeatedly rely on the following consequence of the two properties 
above: Fluctuations of a under ji^ are stochastically controlled on all possible scales by 
arrivals of ^ under the afore mentioned Poisson measure P^. 



3. Main Results 
A slight generalization of our Hamiltonian of interest is given by 

N \ N 

^L^l'^UU^'^t^+hal^^] (3-1) 

where V{-) is some polynomial function from [—1, 1] to R. 

As is apparent from ( |1.7| , an analysis of the phase diagram of the corresponding gen- 
eralized CW model in transverse field boils down to an investigation of asymptotic prop- 
erties of the sequence of measures 





expjN j^VimN{t))dt^ ; da 




'ex^[N j^V{m^{t))dt] 





where 

\ 

This problem may be addressed with the theory of large deviations. 
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Theorem 3.1 For every h > and A > the law of the average m^{-) under the product 
measures is exponentially tight on L2(S^), and, furthermore it satisfies on the latter space 

a LD principle with a good rate function I^j^. 

This follows from fairly standard methods and will be proved in the next section. 
Consequently, the measures ( |3.2| are also exponentially tight on L2(S^) and they satisfy 

a LD principle with a good rate function max^ (5^ j^{m) — (3^ where. 



Khi^) = jy{ni{t))dt - l{^{m). 



In particular, in the zero magnetic field case of = 0, measures (3.2 1 are asymptotically 
concentrated around solutions of 

sup I [^P{m{t))dt - Jf(m)| = sup(&{{m). (3.3) 

Below we shall use (•, ^o denote the scalar product in ^2(6^) and || • ||/5 to denote the 
corresponding norm. Then, since we formulate our large deviation principle in L2(S^), 

= sup Ug,m)f,-4,{g)} (3.4) 

g 

where 

Khis) =log^[^^'''^''']- (3-5) 
For the solution of the variational problem and, in the computation of the corresponding 
transitional curve, one must proceed on a case by case basis. Let us specialize to the 
quadratic problem. 

In order to study phase coexistence in this model, we must consider the case of zero 
magnetic field in z-direction, h = 0. Let us define 

f(A,/3) = ^Var?[((7,l)^] (3.6) 

where Var^ is the variance under the one-circle spin measure fi^. A straightforward 
calculation (which we present in Section |6]l shows that 

f(A,/3) = itanh(A/3). (3.7) 
For the statement of our main theorem, let 

S4(A,i6) = -^A$(/2-l) and S4 (A,/3) = s° (A,/5) . (3.8) 
Our analysis below will show that there exists C = C( A, /3)>0 so that 



d/z3 



ih-l)< -Cho. (3.9) 

h=h„ 
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for ho close to zero. Thus, in view of the spin-flip symmetry and continuity of semi- 
invariants S4 (A, j6) in /z, S4 ( A, j6) is strictly positive. 
The results of our analysis are summarized as follows: 

Theorem 3.2 Let V{x) = The variational problem ( |3.3| l at zero magnetic field h = has 
constant maximizers ±m*{A, fi)-l, where the spontaneous z-magnetization m* satisfies: 

(1) //f(A,|6) < I, then m* =0. 



(2) if f (A, j6) > 1, then m* > 0, and, consequently there are two distinct solutions to (3.3 1 



Furthermore, away from the critical curve the solutions ±m*-l are stable in the following sense: 
There exists c = c{A,(i) > and a strictly convex symmetric function U with a U{r) ~ r log r 
growth at infinity, such that 

iS{{±m*-l) - &{{m) > cmax|||m±m*-l||^,^^LJ(m'(f))df|, (3.10) 

where, 

||m±m*-l||^ = min I ||m — m*-l||^ , ||m -|- m*-l||^| . 

As a result, the variational problem ( |3.3| l is also stable in the supremum norm || ■ ||sup- Namely, 
there exists a constant Cgup = Csup(A,j6) < 00, such that 

«^'^"--"-«AW>exp{- ||„J;Ti||..,j - 

Finally we have the following expression for the decay ofm* > in the super-critical region near 
the critical curve: 

where the implicit constants depend on f) and A but are bounded below in compact regions of the 
parameter space. 



Remark 1. We remark that the second term of ( 3.10| l is particularly important in the 



super-critical regime (f ( A, /3) > 1) since it rules out trajectories of ( ■ ) with rapid tran- 
sitions between the optimal values ±m*. 

Remark 2. We remark that there are alternate approaches to quantum Ising systems in 
which the random cluster representation (and percolation therein) play a more central 
role, e.g. as in [19] [18] and [[T6|. While these ideas will provide a backdrop for our 
derivation, at the core we shall adhere to a more traditional thermodynamic/large de- 
viations approach: Various results, e.g stability, are stronger in this context. As a bit of 
foreshadowing to future work, we remark that the spin system formulation is at present 
essential for the extension to finite dimensional systems. On the other hand we would 
like to stress that our results automatically imply that, as it was claimed in |fl9| and later 
on conjectured in [16[; f(A, /3) = 1, is the equation of the critical curve for the q = 2 
quantum FK-percolation on complete graphs. 
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In Section|4]we shall prove Theorem 3.1 Then the remainder of this paper develops 



our proof of Theorem 3.2 



4. Large deviations under product mesures 

In this section we study asymptotic properties of the average wiAf(-) und er the se- 
quence of product mesures ^fi^fi- Let us begin with the proof of Theorem : 



3.1 



Proof of Theorem 3.1 Let us start with the exponential tightness on L2(S^). Recall the fol- 
lowing realization of compact subsets of L2(S^): Let a : [0,fi/2] K+ be a continuous 
non-decreasing function with a[0) = 0. Then, 

/C„ = UgL2(S^) : ||^||^<land||^(s + -)-^(-)ll/5<«(s)VsG [0,/3/2]} (4.1) 

is compact. Now, in view of the Edwards-Sokal representation, 

1 N 

||mN(s + -)-mN(-)ll/i < 2s-X]^,(S/5). 

(=1 

The latter expression is increasing in ^ = (^i, . . . , ^n), hence it is exponentially bounded 

in view of stochastic domination by product Poisson measures (^iP^. 

Large deviations with good rate function 7^ h defined in ( |3.4| | is a standard consequence 

0. ' □ 

In the sequel we shall need the following stability estimate which implies that product 
mesures (^f^^i, are sharply concentrated around constant functions. 

Theorem 4.1 Let A,h > be fixed, and let rj be an intensity of the Poisson process of arrivals 
on such that ( |2.3| l holds. Then, 

, foo, if m is not absolutely continuous 

~ \loU,^im'{i))d^' otherwise, ' ' 

where is an even smooth strictly convex function with a superlinear growth at infinity. 

Remark 3. The function is defined in ( |4.6| l below. Note that it is negative in a neigh- 
bourhood of the origin. What is important, however, is a super-linear growth at infinity 
which, as we shall see in the concluding Section [7| controls rapid oscillations for the 
original variational problem ( |3.3| . 



Proof. Let TZ = {0 < ti < t2 < ■ ■ ■ < t„ = to < fi} he a partition of S^. Consider an n- 
dimensional random vector 

Zn = [mN{ti) -mN{to),mN{t2) -mj^{ti),...,mM{tn) -m]^{tn-i)]- 
Evidently, satisfies a LD principle with rate function I^, 

l^{zi,...,Zn) = max <J^giZi-A'^{gi,...,g„)\, (4.3) 

Slf-'Sii I ; I 
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where 



(4.4) 



By contraction 1^ 

> I^[m{h)-m{to),...,mitn)-m{tn_^)]. (4.5) 

It, therefore, remains to derive an appropriate lower bound on J^. Alternatively, we 
may seek for an upper bound on A^. Notice that in view of the Edwards-Sokal repre- 
sentation. 



1 + 1 



mi-iM)>o} 



At this stage we are delighted since the latter expression is monotone in ^. Hence, by 
virtue of (pfsb. 



Consequently, 



XigMti)-''it,-i)} 



\ti — ti-l\ 



where H{g) = {e^g + e^^c?). By duality 

I^{Zi,...,Zn) > J2\ti - ti_i\U^ 
i 

with, 

LJ,;(z) = ?/H* , 

with H* being the Legendre transform of H. Thus, is an even smooth strictly convex 
function. Furthermore Utj{m) ~ \m \ log |m| at infinity. Note, however, that U,j is nega- 
tive in a neighbourhood of the origin, in particular 11^ (0) = —2rj. Going back to \4l.5\ we 

obtain that the rate function ,. satisfies 



(4.6) 



i 

for every partition 7?. of S^. Taking supremum over TZ-s we recover ( |4.2 



m{ti) -m{ti_i) 



□ 



This is a somewhat formal argument which, however, could be easily turned into a rigorous one using a 
standard molifying procedure 
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5. Duality and Reduction to One Dimensional Variational Problem 

Recall that we are restricting attention to the case of quadratic interaction 'P(x) = 
II. In particular, the functional 

= \ ^^^m\t)At = ^\\m\\j 

is convex and Gateaux differentiable on L2{Sp,dt). Since, as we already know, has 
compact level sets, the supremum in ( |3.3| l is attained. Let m be a solution, 

Y(m) - I^{m) = max |Y(m) - (m)| 

Then, 

I^{m) — I^{m) > Y(m)— Y(m) > {m — m,m)^, 

for every m G L2{Sp). It follows that is sub-differential at m and, moreover, m E 
dl^{m). But that means A^(m) + I^{m) = Hence, for every m and h, 

A$(m)-^||m||^ > = - (m)} - - Y(m)} . 

We conclude: If m is a solution of ( |3.3| l, then 

Ki^)-^\\m\\j = max|A5(?z)-l||/z||2| (5.1) 

Since is Gateaux differentiable and, in particular, it is ever}rwhere sub-differentiable 
the reverse conclusion is equally true: if m is a solution to ( |5.1[ ) then it is a solution of 
( |3.3| l. Note that far reaching generalizations of the above arguments may be found in 
f27J. An application for general interactions V will be attended elsewhere. Here we 
continue to stick to the quadratic case. Also, since some of the computations below will 
be done in a greater generality, we shall casually rely on the fact that m is a solution to 
( |3.3| l if and only if it is a solution of the dual problem ( |5.1| |. The latter happens to be more 
susceptible to analysis. 

Our next step is to assemble facts which reduce ( |5.1| to a study of constant fields h = c -1 
for c G R. In what follows we rely only on the fact that the polynomial function V has 
a super-linear growth at infinity. In particular, classical Legendre transform V* is a well 
defined function on R. Let integral functional Y : L2(S^) — > R be defined by 

Y(m) = / P(mt)dt. (5.2) 

JO 

and for /z G L2(S^,di) let 

Y^h) = sup{{h,m)^-Y{m)} = / V*{ht)dt (5.3) 
be its Fenchel transform. 
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Lemma 5.1 For any polynomial 'P{x) of a super-linear growth at infinity , let Y and Y* he 
defined as above. Then we have 



sup 



A{{h)-^*{h)] 



sup^ A^(c-l) 



(5.4) 



In particular this lemma implies that the constant function c • 1 optimizes the left hand 
side if and only if c optimizes the right hand side. 

we need the following observation: For any h G L2(S^), let hi 



To prove Lemma 



5.1 



and hr be the pair functions obtained by reflecting h about and That is let 

(h{t) ifie[o4] 

ifiG[f,/3] 



(5.5) 



and hr be defined analogously. 

Lemma 5.2 Let h G ^2(8^) and hi, hy he as ahove. Then we have 

^{{h)<\[A{{hi)+A{{hr) . 

Consequently, for every h e L2(S^), 

< J A{{h{t).l)dt. 



(5.6) 



(5.7) 



We note that Lemma 5.2 was originally proved in a somewhat more general (but com- 
pletely analogous) context in [9J. 

Proof. Inequality \5.6) follows from reflection positivity of classical ferromagnetic Ising 

systems: Recall that is the continuum limit of a sequence of discrete Ising models 
with asymptotically singular interactions. Let M G N be fixed and suppose that the field 
h E V, that is h is piecewise constant on the dyadic intervals [i^f!i, ^^^). Expressing the 
discrete Ising model on a circle via transfer matrices and using the generalized Cauchy- 
Schwartz inequality [12J, one may check equation ( |5.6| l directly whenever the field h G 
D. Passing to limits with respect to the discrete Ising models, we immediately have 
( |5.6| l for any dyadic piecewise constant function h. Standard approximation arguments 
extend the bound to the general case. Finally, \5.7) follows from repeated applications of 

□ 



(5.8) 

□ 



Let us go back to the proof of Lemma 5.1 By ( |5.3[ ) and ( |5.7| |, 

A{{h)-^*{h) < {^^A{{h{t).l)-V*{h{t))^dt 



(|5.4ll follows. 



12 CHAYES ET. AL. 

6. Analysis and stability of the one-dimensional problem 
Following ( |5.4| we need to analyze the following one-dimensional problem. 



sup q{{c) = sup \ A?(c-1) 
Evidently, |a5(c-1)|<j6|c|. Hence, 



(6.1) 



(6.2) 



so that solutions to ( |6.1| always exist. Any such solution c is a critical point satisfying 

(6.3) 



By the FKG and GHS inequalities (cf. the discussion at the beginning of Section 2) the 
function ^M(c) = ^Var^^ (((7,1)^) is positive and decreasing in c. In fact, as we shall 
check below it is strictly decreasing on [0, oo). Therefore Equation ( |6.3| l has nontrivial 
solutions if and only if 



dc 



M(c) 



c=0 



> 1. 



(6.4) 



As {{a, 1)^)'^ = by spin flip symmetry, Var^ \^{(r, 1)^J = f/^ ( [a, 1)^ ) . 
We are reduced to computing 



(6.5) 



where we have used rotational invariance in the equality. In view of the Edwards-Sokal 
representation, the ±1 symmetry implies 



(6.6) 



f{ [im] 

where, as in ( |2.2| before, #(^) denotes the number of connected components determined 
by the underlying Poisson process and {0 denotes the event the and t are in the 
same connected component of the complement of the arrival points. 

There are two cases to consider: either there is an arrival in [t, j5] or there is not. Taking 
this under consideration, we have 



0^i;2*(?) 



,A(/3-2t) 



+ e 



-A(/J-2f) 



A simple computation shows 



2#(S)" 



-A/5 ^ gA/i 



(6.7) 
(6.8) 
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Upon integrating, we find 



Var;((,7,l)s) = -lanh(A/!). (6.9) 



Conclusions (1) and (2) of Theorem |3 . 2| folio w as soon as we show that 

1 

is indeed strictly decreasing on [0,oo). To this end, note that. 



dc 

where. 



lo L L ^A,c(^S'Odrdsdt, (6.10) 



- <^W^(0>l,(^(s)>t + 2(^(r)>^^^(^(s)>5^^(c.(0>,^^, (6.11) 

the third Ursell function for In view of the above mentioned identification of ji^^ 
as a limit of one-dimensional ferromagnetic Ising models, we state the following result 
without proof: For any c > 0, 

U{^^{r,s,t) < forallr,s,f G S^. (6.12) 

In order to derive a strict bound we shall employ a similar approximation, which will al- 
low us to transfer the random current representation used by Aizenman [1 J and Aizenman- 
Fernandez [4] to this continuous time setting. We obtain the following lemma as a result. 
While the proof is not difficult, it is tedious and so we relegate it to the appendix. 

Lemma 6.1 [Random Current Bound] For any -positive constant field c and any triplet of points 
r,s,t G S^, we have the hound 

Remark 4. We note further that one can show that 

<,(r,s,0 > - C(A,/3) (A/3)2 {a,) (6.14) 

for some C > 0, bounded away from zero on compact regions of parameter space. This 
will not play any role below and so we do not pursue the bound here. 

Thereby, the phase diagrame ( |6.4| | is justified. Let us turn to the stability issue: As 
follows from ( 6.10| and ( 6.13 1, there exists X = Xi^r ^) r such that for any c > 0, 



Ul{r,s,t) < -{^o)le-(''^^^)^^ [{s - rf + {t - sf + {fi + r - t)^] . (6.13) 
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As a result, for any < ci < C2, 



A,C2 



(6.15) 



Assume that ( |6.4[ | holds, and let c* = m* > be the positive maximizer of 0^. Then, of 
course, = ^0^(c*). Let us derive an upper bound on 



For c > c*, 

-I 

c Jo 



dc 

-.1 

c* Jo 



T — 

L C* 



:* Jo I 



dr = 1 + 



T- 



4{T)}dT. 



By ( 6.15| the second term is bounded above by 



X 



rdr 



-Xic-c*)- 



c* 



{croY^TdT < 



X{c 



(c-c*), (6.16) 



where at the last step we have relied on (c^o)^ t — ^ every t < c* 
In a completely similar fashion we deduce that, 

>r<(x)dx-i>5(^(, 

c Jo 3 
whenever, c < c*. We have proved: 



(6.17) 



Lemma 6.2 Assume ( |6.4| |. T/ie f?ze positive maximizer m* = c* > of Q)^ is stable in the 
following sense: For every c G [0, oo). 



(6.18) 



7. Conclusion of the Proof of Theorem |3]2j Claims ( |3.10[ ) and ( |3.12[ ) 

Let us return first to the infinite dimensional variational problem ( |5.1| |. The FKG prop- 
erty of ji^ (inherited in a limit of ferromagnetic Ising models) evidently implies that 
A$(/2) < Aj(|/2|).Sinceby g. 



Km 



< 



PJO 
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we readily recover from (6.18 1 the following L2-stability bound: 

{a$(..i)-M-:)^}-{aJ(|.|)-^IHI|,} 



- / d{\h{t)\-m*)dt ^ D{h). 
p Jo 



> 



In order to transfer this bound to the original variational problem ( |3.3| l we shall again 
rely on duality considerations of Section |5] Recall that 

A$(m*-1) - /3(m*)V2 = /3(m*)V2 - (7.2) 



Therefore, by ( |7.1| | 

/3K)V2-7jK-l) > {{h,m)^-\\h\\j/2}-I^^{m) + D{h). 
Taking h = m we arrive at 

©^(m*-l) -©^(m) > D(m). (7.3) 

Obviously, there exists Ci = Ci (A, j6), such that 

D{m) > ci\\m±m* ■ 1||^. (7.4) 

Thereby, the first term in ( 3.10| l is recovered. In order to recover the second term in (3.10 1, 
set in = m*-l and just compute: 

1 1 

+ I^{ni) — I^{fn) — {m,m — fn)^ 

= --||m -m||^ + lj^,„*(m). 
Repeating the same computation with — m, we infer: 

<d{{m)-(5{{m) > -^-\\m±m\\l + min{ll^^{m)jl_,^4m)] 

1 fl^ 
> --\\m±m\\j + U^(^^^„.){m'{t))dt, 

as follows from (4.2 1. At this point we recall ( |7.4| |, and the full statement of ( 3.10| l becomes 
an easy exercise (with U = U,j). 

Let now m be such that ||?« ± m||sup = 2^. Let us consider only the case when 
maxf = m* +2S. lfmintm{t) > .5, then \\m ±f«||^ > S^fi- Otherwise, without loss of 
generality, we may assume that m (0) = m* +23 and let r = min {t > : m{t) = m* + 3}. 
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Obviously, in such circumstances, 

\\m±m\\j>rS^ whereas U{m' {t))dt > j'^U{m' {t))dt -2rj^. 
However, by convexity. 



U{m'{t))dt > rU{S/r) ~ S\og{l/r) 



for r <C S. Therefore, the J^^ LJ(m'(i))dMerm in max | ||m ± Jp^ L[(m'(i))di| be- 
comes dominant for r ~ e^^^^. Hence ( |3.11| . 

Finally, consider the claim (3.12|. As we have seen in Section|5| m*!-l solves the dual 
variational problem \5.1\ . In particular, 

1 fi 



m 



^ ^ \,m* 



Via the Taylor expansion we may write 



m 



^Var^,,((^,l) 



*^3 



S4 (A,j6) K ) 
6/3 



+ 0(K)4 



(7.5) 



(7.6) 



noting that odd derivatives of A^ vanish at zero by spin flip symmetry. As we have 
already mentioned. Lemma 6.1 implies that S4(A, j6) > 0. Claim ( |3.12 l follows easily. 

□ 



8. Appendix 



Proof of Lemma 6. 1 We begin by recalling the random current representation of the third 
Ursell function for a discrete Ising model, see [4]. 

Following the notation of HI, consider a general finite graph Ising model on a vertex 
set V (also called sites), elements of which we denote with i, and with bonds b in some 
fixed subset 6 C V x V. Let us denote the coupling constants by and local fields by 
hj. It is useful to interpret the local fields /z, as the coupling constants between the site i 
and a 'ghost' site q. We augment the graph V, £ by including g along with an edge set 
connecting q to each vertex in V. Let us denote the augmented edge set by £'. 

Let n = {ni,)},(z£i denote a sequence of integer valued 'fluxes' attached to the bonds 
of the underlying graph. We say that that a site z S V is a source if J^i^h rii, is odd and 
denote the collection of sources other than the ghost site by dn. 

We say that x ^ y \{ there exists a path of non-zero fluxes connecting x to y using 
bonds in £ (i.e. not ghost site bonds). Moreover, x h means that whenever x ^ y, 
riy^g = 0. For any subset of bonds B C £, let 

(■>H (8-1) 
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denote the Gibbs state for the Ising model with the constants {}b)b(^B *° ^- 



W(n) 



Ju"'' 



and 



(8.2) 



(8.3) 



z = ^ Win). 

dn=0 

Suppose r,s,t G V are fixed sites in the Ising model under consideration. To avoid 
excessive provisos hereafter we shall assume that all three points are distinct. For a 
collection of fluxes n, let us say that b S C„(r) if z <-> r for some i G b. Then we have (the 
reader should consult for a derivation) 




W(ni) W(n2). 



h X 



Qj+„2('-) 



+ {s^t} (8.4) 



where {s <4> represents the first term with the roles of s and t interchanged. Here 
and are independent copies of fluxes. 

Clearly, the weights W(n) are proportional to the probability that a family of inde- 
pendent Poisson processes indexed by (generalized) bonds take a collection of values 
determined by n. We need to differentiate between the processes associated to the bonds 
of the graph and the bonds with the ghost site g. Specifically, let {Afh, -MiybeSj^v denote 
a collection of independent Poisson processes with respective parameters 
{Jbfhi}b(^£j(^v srid let dP denote the joint probability measure associated to these pro- 
cesses. 

It is well known that if h > 0, spin correlations are increasing with respect to cou- 
pling strengths, so each summand on the righthand side of ( |8.4| must be non-positive. 
Therefore, neglecting summands we obtain 



Qj+„2(»-) 



+ {s^t} (8.5) 



W{n,)W{n2) 

where we drop analogous terms in the expression {s <4> f }. 

Let us concentrate on the first term on the righthand side of the inequality as the 
second may be treated by symmetric considerations. We note two things. First, on the 
set [ni + : r ^ n {n^ -^n^'.r ^ we have 

l^^^^^ = (8.6) 



18 CHAYES ET. AL. 

by spin flip symmetry. Second, we have 

{dni = {r,s},dn2 = 0} fl {n^ + n2 : r ^ h] = 

{n^:r^s}r\ {dn^ = {r,s},dn2 = 0} n {n^ + n2 : r ^ h}. (8.7) 

In terms of the Poisson processes, let Er^s,t denote the event determined by the require- 
ments of the first term on the righthand side of (8.5 1. After the reduction made above, 
( |8.5| l may be expressed as 

<(r,s,f) < + {s^t} (8.8) 

(P(ni = 0)) 

Let us now specialize to Ising models on the circle with 2^ vertices. We take coupling 
and field strengths given by 



For the readers convenience we note that 



1 _ M 

P(A4 = odd) oc tanh(/,) = 1^ 

P(A4 = even) oc 1 
P(A/-fc = 0) occosh(/fc)-i = yp + 0(2-~) 



and 



F{Mi = odd) = ^e-l^ +0{2-^^) 

F{Mi = even) = e ^ + 0{2-™) 

_ hf, 

V{Mi = 0) = e 2~ 
where the implicit constant of proportionality is 

(l + o(l))e-"''2~. (8.9) 

We compute a lower bound for the probability determined by the numerator in ( |8.8| |. 
Without loss of generality we may assume the orientation 

< r < s < i < j6. (8.10) 

Continuity arguments imply that it is sufficient for us prove the bound ( 6.13| l assuming 
each point is of the form /'//2* for some /, A; G N fixed. 

Let c,d e [0,j6] be fixed. Let [c,d\ denote the integers i in [0,2^^] such that ^ < / < 

If c > d this is understood to be with respect to periodic boundary conditions. For 
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any bond & let G [c, d\ denote the case that both vertices of h are in [c, d\ . Finally, let 

M^cdj = {Mf^ + Mf^ > for some i G [c,d\ }. (8.11) 

Based on topological considerations there are three cases: £r,s,t H M^r,sJ/ ^^/-,s,t H M^, fj 
and Er^s,t H M|^t^j . The requirements of Er,s,t imply these three possibilities lead to disjoint 
subevents (but not a partition) of Er,s,t- 

The cases carry a minimal but tedious amount of computation, so we shall present 
only one of them. 

Let us consider Er,s,t n M^f ,.j . For any x<r<t<y in the dyadic lattice {^} ,\et 

A{x,y) = {A/;*^' +Arf^ > Oforallfc e [x,r\ U [t,y\}n 

{MI'^ + A/'f ' = 0, if & = {(--1)2%, -2~/^} or {y2%, (y+i)2%}}. (8.12) 

The parity of all b-bonds in [x,y\ is completely determined on the event Er^s,t H 
M[^f n A(x,i/). Since the fluxes associated to g-bonds must be zero, independence of 
the various Poisson processes allows us to compute: 

P {Er,s,t n M^rAi n M^^y)) ^ e-2^l-'-|-2/.|i/-x| + 0(2-™). (8.13) 

A worst case lower bound for the exponential term is e-(2A+2;i)/s^ gumniing over or- 
dered pairs {x,y} and taking into account ( |8.9| , 

P {E,.,s,t n ML,,fj ) > e-(4A+2/<)/3 1 [A (/3 + r - f)]' + o(l). (8.14) 
Analogous estimates in the other cases lead to 

P {Er,s,t) > ^e-(4A+2/^)/i [(s _ rf + (i - s)2 + + , _ i)2] + 0(1). (8.15) 
As P(9n = 0) < 1, the bound ( [6l3| follows. 

□ 
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